Unruh gamma radiation at RHIC? 
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Varying the proposition that acceleration itself would simulate a thermal environment, we investi- 
gate the semiclassical photon radiation as a possible telemetric thermometer of accelerated charges. 
Based on the classical Jackson formula we obtain the equivalent photon intensity spectrum stemming 
from a constantly accelerated charge and demonstrate its resemblances to a thermal distribution for 
high transverse momenta. The inverse transverse slope differs from the famous Unruh temperature: 
it is larger by a factor of tt. We compare the resulting direct photon spectrum with experimental 
data for Au-Au collisions at RHIC and speculate about further, analytically solvable acceleration 
histories. 
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THERMAL LOOKING CLASSICAL RADIATION 
FROM ACCELERATING SOURCE 



It has been a puzzle to understand why thermal mod- 
els supposing global equilibrium, and hydrodynamics as- 
suming local equilibrium are so successful in describing 
single particle observables in high energy heavy ion and 
even in elementary particle reactions. Neither is it easy 
to understand intuitively how accelerating frames pro- 
duce " thermal vacua" , with quanta distributed as if they 
were connected to a heat bath - a proposition first for- 
mulated by Unruh 1] and soon analyzed in relation to 
black hole thermodynamics by Bekenstein ^ and Hawk- 
ing ji, JJ . Spectra of photons from a black hole was ob- 
tained by Page These early analyses in the seventies 
were followed by numerous attempts to harvest the mer- 
its of this idea in relation with possible ways to quantum 
gravity [6j . In the present letter we calculate measurable 
consequences of this idea in photon spectra observed in 
high energy heavy ion collisions. The observation that 
transverse momentum spectra of direct gammas, mea- 
sured perpendicular to the beam line in collider experi- 
ments, are well fitted by exponential distributions 7], are 
suggestive of " thermal" heat bath sources in high energy 
nuclear collisions. Theoretical calculations assuming the 
applicability of thermodynamic concept are phenomeno- 
logically successful in describing heavy ion experimental 
data [8l4ld| . Statistical predictions based on the assump- 
tion of a thermal source for also pp collisions at the LHC 
have recently been made by Becattini et al. 
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Alternatively, the proposition that the thermal spectra 
'mp + p and in heavy ion A + A collisions were due to ac- 
celerating frame concepts akin to the Unruh temperature 
of field theoretic vacua was made by Satz and Kharzeev 



However, a dynamical calculation demonstrat- 
ing how thermal looking spectra emerge from accelerat- 
ing systems was not provided. 

In this letter we demonstrate how the classical elec- 
tromagnetic field of a constantly accelerating charge de- 
velops an approximate thermal distribution in the trans- 
verse photon spectrum while also exhibiting a uniform 
rapidity dependence usually interpreted in terms of lon- 
gitudinally boost invariant Bjorken hydrodynamics. This 
calculation is entirely classical with the Planck constant 
entering solely by converting the intensity distribution in 
terms of the number distribution of equivalent photons 
as each carrying the energy huj. Beyond this we do not 
need to refer to field theory, quantization procedure or 
the vacuum state at all. The trajectory of pointlike clas- 
sical charges and the radiated photon field are calculated 
in the framework of special relativity, without any actual 
reference to a general frame and related Killing vectors. 

We consider a simple point charge undergoing constant 
acceleration. We calculate the classical bremsstrahlung 
distribution using the text book[23| formula for the radi- 
ation field amplitude: 



A = K 



dt 



n X {n X (3) 
l-n/3 



dt 



(1) 



with fC^ = e^/ 
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and retarded phase ip = uj {t — nf/c) . 



Here the trajectory of the moving charge is given by r(t) 
and (3 = V / c = dr/cdt. The relativistic motion of the 
charge can also be described in terms of a proper time r. 



This defines a Lorentz factor 7 l/yl — /^^ = dt/dr. 
This way the derivative of the retarded phase with re- 
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spect to the proper time is given as 
^ = W7(l - n(3). 



(2) 



Without restricting the generahty the velocity vector can 
be taken to point into the first spatial direction, /3 — 
(w/c, 0, 0), with the magnitude v/c = th^, ^ being the 
rapidity variable in the direction of the motion. The 
original Unruh formula arises by a Fourier analysis of 
the phase observed in the particular detecting direction 
of zero angle. 

In our calculations we use the vector 



n X (n X /3) 
1 - n/3 



V sin 



c — V cos f 



(3) 



with e — dfi/dO unit vector orthogonal to n. The Doppler 
factor is given as 



7(1 — (in) — ch^ — sh^ cos ( 



(4) 



Next we define thry = cos 9. Due to the extreme rela- 
tivistic kinematics 77 = lnctg(0/2) is the rapidity of the 
photon. By the virtue of this, one arrives at the following 
Doppler factor: 



7(l~/3n) 



ch(g-7?) 
ch77 



(5) 



Since this is proportional to the derivative of the retarded 
phase with respect to the source's rapidity. 



dip g dtp ch(^ — 77) 
dr c d^ chrj ' 



(6) 



by integration with constant acceleration, we arrive at 
Loc sh(^ — 77) 



g chr] 



(7) 



Considering that lu = ck±/ siiiO ~ cfcj^chyy and I — (? j g 
is a length characteristic to the magnitude of accelera- 
tion, the above result can be written in a simpler form; 



[p — ik± sh{£_ — if). 
The vector u can be expressed as 

sh$ 



/3sin6' ^ 
u ~ e 



1 - /? cos 6* " ch(C - ri) 
its derivative is given by 

du chq 



(8) 
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The Lorentz covariant radiation intensity distribution is 
evaluated here in terms of transverse wavenumber and 
rapidity variables. The correspondence to the familiar 



covariant notation is established by considering the fol- 
lowing four-vector definitions: 

(7,7^^) = (cK,she,0,0) (11) 

for the four- velocity of the radiating charge, 

k = [uj,ujn) = /s^(ch?7, sh?;, cos 7/;, sin?/;) (12) 

for the four-momentum of the photon and 

1 dk 

ei = — — = (0, — sin0, cos^cost/i, cos^sinT/!), 
LO oO 



£2 = 



1 dk 
uj sin 6 dip 



— (0, 0, — sin -0, cost/'), 



(13) 



for the polarization vectors. Using these covariant nota- 
tions the familiar difference for the amplitude is replaced 
by a continuous integral over rapidity 



e • u 



k-' 



d 



I d£^ \k ■ u 



dt 



(14) 



The Jackson formula ([T]) is more in taking into account a 
continuously changing phase along the path of the mov- 
ing charge, arriving at 



d^ 



Finally the spectrum is given as 
dN 



kj^dk±dri dtp 



2\A\' 



(15) 
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In any case the description of the linear motion of the 
source of the radiation is given by the proper acceleration, 
d^/dr = g{£,)- Accordingly the phase has to be obtained 
from the integration of the relation 



dip d ^ 
dr dr 



k ■ u = fcj^ch(^ — 77). 



(17) 



The general trajectory with constant relativistic accel- 
eration, the Unruh trajectory, is given parametrically as 



i = to + -(sh(e + eo)-sh$o), 
9 

2^3 = 2^03 H (ch(^ + ^0) - ch^o) • 

9 



(18) 



with (3 = v/c ~ th(^ -f ^0) and ^ = 9t/c. The phase 
space factor of the radiation can be expressed by using 
the fc_L — Ljsm9/c transverse momentum and the rapid- 
ity 77: One obtains diudn — c{dk j^drf) / chrj . This leads 
to the following rewriting of the differential spectrum for 
the photon number: 



dN 



k± dk± di] dtp hk\ ch r] 



A 



(19) 
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By determining the amplitude of the radiation it is note- 
worthy that the derivative of the vector u with respect to 
time and the integration with respect to the same time 
can be replaced by any variable. The amplitude is given 
by using the above expressions for ip and u as functions 
of^: 



+ 00 



A^Ke 



chr; 



From this one obtains the following analytic result 

A = 2Ke— Ki {Iki^) 
9 



(20) 
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with the Bessel K-function. The spectrum of equivalent 
photons is given by 



dN 



^OLEM 



k±dk±dr] dip 



(22) 



with uem — e'^/A'Khc. 

For heavy ion collisions basically two limiting approx- 
imations can be considered: the coherent square of sum 
of amplitudes and the incoherent sum over squared el- 
ementary amplitudes. For k-r > I GeV and i? ~ 7 fm 
nuclear radii, coherence is negligible and the incoherent 
sum leads to an effective charge factor 



F = 2Z{b) 



(23) 



with Z{b) being the proton participant number of each 
nucleus colliding at impact parameter b. Z{0) = 79 most 
central & = for Au beams, while (Z) « 55 for the 
— 20% centrality trigger bin considered in the PHENIX 
experiment. 



THERMAL INTERPRETATION: NOT EXACTLY 
UNRUH 



For low k± the spectrum (|22p resembles the conformal 
l/k\ limit independent of I. For high fc^ on the other 
hand it develops an approximate exponential tail resem- 
bling thermal fireball models: 



dN 



k±dk±dri dip 



F 2aEM 



e 



-2lk, 



(24) 



The equivalent temperature is seen from the exponent 
and can be interpreted if we wish as tt times the Unruh- 
temperature, Tj/: 



Tic 



2kBt 2kBc" 



ttTtj 



(25) 



Note however that this photon spectrum is uniform in 
rapidity resembling Bjorken hydrodynamic scenario. The 
uniform rapidity distribution is nevertheless truncated at 



large 77 if we consider a finite range of time for which g 
is constant. 

Although the photon spectrum obtained above does 
contain an exponential tail, and defines an effective trans- 
verse temperature, it is clear from this specific dynami- 
cal model that no thermally equilibrated source was ever 
involved. It only "appears as if" a locally equilibrated 
Bjorken hydrodynamic fluid sourced the photons. While 
the classical radiation spectrum of photons from a con- 
stantly accelerated point charge shows an exponential 
tail, and by the virtue of this defines a temperature re- 
lated to the acceleration, g, the details of the resulting 
photon rapidity and k^_ distribution are not identical to 
that assumed with thermal models. 

Instead of a thermal Bessel Kl function, its square 
occurs, and instead of the Unruh temperature its tt times 
enlarged value governs the tail. It is therefore important 
to understand better the classical EM result and the role 
of the K-Bessel function. The latter emerges not only 
from the EM Fourier integral above, but also from the 
relativistic statistical Gibbs distribution. Also the angle- 
and rapidity variables, and integrals over them, occur in 
relativistic calculations with an assumed collective fluid 
velocity fleld. It is possible to reinterpret eq.(22) in terms 
of a Jiittner-type representation by using the identity 



+00 



dC, K2 {z chC ) = Kl 



(26) 



Utilizing this relation we get 



dN 



k± dk± dr] dip 



+00 

^ioE^ Jd(K2i2ik^cHC-r^)\ 



This shows that ordinary bremsstrahlung for a uniform 
accelerating charge can indeed be thought of a superpo- 
sition of fluid cells distributed with uniformly distributed 
collective fluid rapidity, referred to longitudinal flow, 
evaluated at a flx freeze-out proper time when the local 
temperature has cooled to eq.(25). It is rather remark- 
able how two completely different dynamical frameworks, 
classical EM and classical hydrodynamics, can lead to 
closely identical "predictions" or more properly identical 
flts to spectra resembling eqs. (22,27). 

Finally we consider a finite period of uniform accelera- 
tion. In that case the rapidity integration is done between 
finite limits: 



dN 



k^dk^drj dip 



= F 



aEM 1 

TT kl 



i£k±w 



dw 



(l + u;2)3/2 



(28) 

with wi^2 — sh(fi^2 + ^0 — v)- Numerically a rapidity 
interval from —3 to -1-3 suffices to approach the idealized 
\wi\ — )■ Qo limit. Thus at high energies at RHIC and 
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LHC where the total rapidity interval exceeds 10, the 
end elfects can be safely neglected at central rapidities 
|,7|<^1 

In the low energy opposite limit, on the other hand, 
the integral can be estimated for W2 — Wi small, and the 
integrand evaluated at Wq — {wi + W2)/2 to get 



dN 



OLEM 1 1 



k±dk± drj dtp 



(29) 



with X = (Ci + ^2)72 + ^0 ^ For a symmetric scenario 
of decelerating trajectory x — V ^nd the above formula 
resembles a Landau hydrodynamic scenario in terms of 
the photon energy: 



d^N 



d^N 



d^k uiduidfl 



F 



O-EM 



1 



TT w^ch ?7 



(30) 
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FIG. 1: (Color online) Acceleration profiles against the ra- 
pidity of a point charge. The photon spectrum is analytically 
calculable for these cases. 



COMPARING VARIOUS ACCELERATION 
PROFILES 

For comparison we consider several simple analytically 
solvable trajectories. One engineered to give a pure expo- 
nential is defined hy g — (7och^(l -I- ch^) and t — th(^/2) 
while v/c — th(^). It leads to a pure exponential spec- 
trum at = 90 degree {-q — Q), 



dN 



k±dk±drj dip 



(31) 



ri=0 



This has, however, an unrealistic acceleration profile, 
having its minimum at ^ = as gmin = '^go while 
£ = c^go. 

A more realistic example is implicitly given by the 
trajectory t = tg(7ru/2) and a; = 1 -I- wtg(7rw/2) -|- 
{2/tt) lncos(7rw/2). In this case the acceleration is mostly 
zero, it emerges up to g^n /2 at ^ = 0. The resulting pho- 
ton spectrum is 



dN 



k±dk±drj dtp 



= F 



Q-EM 1 
TT kl' 



-2£fc 1 



(32) 



Finally we note that there is always a thermal looking 
part of the photon spectrum beyond the transverse mo- 
mentum for which the retarded phase is larger than tt. 
This is achieved for 



AC 



cAt shA^ 



(33) 



if estimating the deceleration by g/cK, A^/At for short 
stopping. It is noteworthy that even for A^ — >■ there 
is a thermal tail, for finite rapidity acceleration histories 
this region is pushed towards lower k± values. 

In Figure [1] we present the proper acceleration as a 
function of the rapidity of the source charge for three 
different, analytically solvable scenarios discussed above. 
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FIG. 2: (Color online) Comparison of the 

Majcwell/Unruh/Bjorken like bremsstrahlung formula 
eq.(22) with RHIC PHENIX data on direct photons in 20% 
central Au-Au collisions for the constant g — 2-kTu scenario 
applying Tu = 135 MeV and the effective charge Z{b) = 55 
(dotted line). The pure exponential fit by the PHENIX group 
fitting the temperature to T = 240 MeV and the amplitude 
to A — 45.2 is indicated by the dashed line. 



The realistic one starts and ends with low acceleration 
magnitude and reaches its maximum at ^ = simulating 
the strongest medium effect in the laboratory frame. 

In Figure [5] we compare transverse gamma spectra 
arising from the simplest (constant acceleration) Unruh 
like scenario bremsstrahlung with Au-Au experimental 
data at RHIC According to the publication cited 
above at low px — hkx the spectrum is dominated by 
a "thermal" contribution, fitted by a pure exponential 
form Aexp(-pT/r) with A = 45.2 and T = 240 MeV. 
For comparison we show the Maxwell/Unruh/Bjorken re- 
sult with the corresponding factor £ = 1/{2ttTu) in the 
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proper units, using the temperature of Ty = 135 MeV. 
Here the magnitude is given once the deceleration fac- 
tor g - or the corresponding Unruh temperature, Tu - is 
fixed. Our primary aim here is not to fit in detail the ex- 
perimental data but simply to point out the degeneracy 
of various theoretical models in terms of bremsstrahlung, 
Unruh , and Bjorken 1-l-lD perfect fluid hydrodynam- 
ics interpretations. The open question is the unique- 
ness of the interpretation of direct photon data in terms 
of purely thermal concepts. How much of the observed 
spectra is due to accelerating classical field effects? Here 
we limited the consideration to pure classical EM ob- 
servables. In the Color Glass Models[26!| large classical 
chromo field effects were also shown to be able to fit 
certain hadronic observables in heavy ion reactions simi- 
lar to thermal/hydrodynamic models. Only through de- 
tailed systematic studies of -y/s, A, and impact parameter 
b scaling in pp, pA, AA of a wide variety of experimen- 
tally observable data can these, in their physical back- 
ground differing, pictures and mechanisms be properly 
sorted out. 



Acknowledgment 

This work has been supported by the Hungarian Na- 
tional Research Fund, OTKA (K68108), by a Hungarian- 
South- African TeT project (TeT-10-1-2011-0061, ZA- 
15/2009), and by the TAMOP 4.2.1./B-09/1/KONV- 
2010-0007 project, co-financed by the European Union 
and European Social Fund. M.G. acknowledges also 
support from Division of Nuclear Science, U.S. Depart- 
ment of Energy, under Contract DE-FG02-93ER-40764. 
T.S.B. and Z.S. thank the support by the Hclmholtz 
International Center for FAIR within the framework of 
the LOEWE program (LandesOffensive zur Entwicklung 
Wissenschaftlich-okonomischer Exzellenz) launched by 
the State of Hesse. Discussions with Prof. I. Racz are 
gratefully acknowledged. 



[2] J. D. Bekenstein, Phys. Rev. D7(1973)2333, 

[3] J. M. Bardeen, B. Carter and S. W. Hawking, Comm. 
Math. Phys. 31(1973)161. 

[4] S. W. Hawking, Comm. Math. Phys.43(1975)199. 

[5] D. N. Page, Phys. Rev. D13(1976)198. 

[6] R. M. Wald: Quantum Field Theory in Curved Space- 
time and Black Hole Thermodynamics, University of 
Chicago Press. 1994 

[7] A. Adare et al. for PHENIX, Phys. Rev. 
LettlO4(2010)132301 

[8] A. Dumitru, U. Katscher, J. A. Maruhn, H. Stocker and 
W. Greiner, Phys. Rev. C51(1995)2166 

[9] J. Cleymans, K. Redlich and D. K. Srivastava, Phys. Rev. 
C55(1997)1431 

[10] T. Hirano, S. Muroya and M. Namiki, Prog. Theor. 

Phys.98(1997)129 
[11] R. Baier, M. Dirks, K. Redhch and D. SchifT, Phys. Rev. 

D56(1997)2548 
[12] P. Aurenche, F. Gehs, R. Kobes and H. Zakaret, Phys. 

Rev. D58(1998)085003 
[13] F. D. Steffen and M. H. Thoma, Phys. Lett. 

B510(2001)98, Erratum-ibid.66O(2008)604 
[14] R. Chatterjee, D. K. Srivastava and U. Heinz, Proc. of 

Eilat 2008, Particles and nuclei (PANIC08), 534 
[15] J. R. Bhatt, H. Mishra and V. Sreekanth, arxiv: 

1005.2756, 2010 
[16] R. Chatterjee, H. Holopainen, T. Renk and K. J. Eskola, 

Phys. Rev. C83(2011)054908 
[17] F. Becattini, P. Casterina, A. Milov and H. Satz, Eur. 

Phys. J. C66(2010)377 
[18] F. Becattini, P. Casterina, A. Milov and H. Satz, J. Phys. 

G38(2011)025002 
[19] D. Kharzeev and H. Satz, Phys. Lett. B334(1994)155 
[20] D. Kharzeev and K. Tuchin, Nucl. Phys. A753(2005)316. 
[21] P. Casterina, D. Kharzeev and H. Satz, Eur. Phys. J. 

C52(2007)187. 

[22] T. Brasovenau, D. Kharzeev and M. Martinez: In Search 
of QCD - Gravity Correspondence, In "The Physics of 
the Quark Gluon Plasma", Eds. S. Sarkar, H. Satz, 
B. Sinha, Springer, 2009. 

[23] D. Jackson: Classical Electrodynamics, Chap. 14, Wiley, 
New York, 1975. 

[24] T. S. Biro and Z. Schram, Acta Physica DebrecinaXLV 
(2011)176. 

[25] A. Adare et al. , Phys. Rev. C81(2010)034911. 

[26] L. McLerran, Acta Phys. Polon. B 41, 2799 (2010). 
F. Gelis, E. lancu, J. Jalilian-Marian, R. Venugopalan, 
Ann. Rev. Nucl. Part. Sci. 60, 463 (2010). 



[1] W. G. Unruh, Phys. Rev. D14(1976)870. 



